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Abstract 

The dipole response of the = 50 nucleus ^*^Zr was studied in photon-scattering experiments at 
the electron linear accelerator ELBE with bremsstrahlung produced at kinetic electron energies 
of 7.9, 9.0, and 13.2 MeV. We identified 189 levels up to an excitation energy of 12.9 MeV. 
Statistical methods were applied to estimate intensities of inelastic transitions and to correct the 
intensities of the ground-state transitions for their branching ratios. In this way we derived the 
photoabsorption cross section up to the neutron-separation energy. This cross section matches well 
the photoabsorption cross section obtained from (7, n) data and thus provides information about 
the extension of the dipole-strength distribution toward energies below the neutron-separation 
energy. An enhancement of El strength has been found in the range of 6 MeV to 11 MeV. 
Calculations within the framework of the quasiparticle-phonon model ascribe this strength to a 
vibration of the excessive neutrons against the N = Z neutron-proton core, giving rise to a pygmy 
dipole resonance. 
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I. INTRODUCTION 



Gamma-ray strength functions, in particular dipole-strength functions, are an important 
ingredient for the analysis of photodisintegration reactions as well as of the inverse reactions 
like neutron capture. These reactions play an important role for specific processes of the 
nucleosynthesis. Moreover, an improved experimental and theoretical description of neutron- 
capture reactions is important for next-generation nuclear technologies. 

The measurement of dipole-strength distributions from low excitation energy up to the 
neutron-separation energy delivers information about magnitude and structure of the low- 
energy tail of the giant dipole resonance (GDR), in particular about possible additional 
vibrational modes like the so-called pygmy dipole resonance (PDR). Dipole-strength dis- 
tributions up to neutron-separation energies have been studied for only few nuclides in 
experiments with monoenergetic photons (see, e.g., Refs. [H, 0, 0|) experiments 
with bremsstrahlung (see, e.g., Ref. jl] and Refs. therein). The bremsstrahlung facility ^ 
at the superconducting electron accelerator ELBE of the research center Dresden-Rossendorf 
opens up the possibility to study the dipole response of stable nuclei with even the highest 
neutron-separation energies in photon-scattering experiments. In the course of a system- 
atic study of dipole-strength distributions for varying neutron and proton numbers in nuclei 
around A = 90 0, @] we have investigated the = 50 nuclide ^°Zr. 

In earlier experiments, the 2+ states at 2186, 3308, 3842, 4120, 4230, 4680, and the 
1*^"-* states at 4580 and 5504 keV were studied with bremsstrahlung produced by 5 MeV 
electrons fol . [lo| . About 20 further states were found in experiments at higher energies up 
to 10 MeV jlll. 1^. For the 15 J = 1 states above 6 MeV negative parity was deduced from 



an experiment with polarized photons [13]. 

In the present study we identified about 190 levels up to 12.9 MeV. We applied statistical 
methods to estimate the intensities of inelastic transitions to low-lying excited levels and 
to correct the intensities of the elastic transitions to the ground state with their branching 
ratios. The dipole-strength distribution deduced from the present experiments is compared 
with predictions of the quasiparticle-phonon model. 



II. EXPERIMENTAL METHODS AND RESULTS 

The nuclide ^°Zr was studied in photon-scattering experiments at the superconducting 
electron accelerator ELBE of the research center Dresden-Rossendorf. Bremsstrahlung was 
produced using electron beams of 7.9, 9.0, and 13.2 MeV kinetic energy. The average currents 
were about 330 /lA in the measurement at 7.9 MeV and about 500 /lA in the measurements at 
higher energy. The electron beams hit radiators consisting of niobium foils with thicknesses 
of 4 fim in the low-energy measurements and of 7 /im in the measurement at 13.2 MeV. A 10 
cm thick aluminum absorber was placed behind the radiator to reduce the low-energy part 
of the bremsstrahlung spectrum (beam hardener). The collimated photon beam impinged 
onto the target with a fiux of about 10^ s~^ in a spot of 38 mm diameter. The target 
was a disk with a diameter of about 20 mm to enable an irradiation with a constant fiux 
density over the target area. The target consisted of 4054.2 mg of ^°Zr02 enriched to 97.7%, 
combined with 339.5 mg of ^^B, enriched to 99.52% and also formed to a disk of 20 mm 
diameter, that was used for the determination of the photon fiux. Scattered photons were 
measured with four high-purity germanium (HPGe) detectors of 100% efficiency relative to 
a Nal detector of 3 in. diameter and 3 in. length. All HPGe detectors were surrounded by 
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FIG. 1: Parts of a spectrum of photons scattered from ^^Zr combined with ^^B during the irradia- 
tion with bremsstrahlung produced by electrons of an energy of Ef'' = 13.2 MeV. This spectrum 
is the sum of the spectra measured with the two detectors placed at 127° relative to the beam. 
The most dominant transitions assigned to ^^Zr are marked with their energies in keV. 

escape-suppression shields made of bismuth germanate scintillation detectors. Two HPGe 
detectors were placed vertically at 90° relative to the photon-beam direction at a distance of 
28 cm from the target. The other two HPGe detectors were positioned in a horizontal plane 
at 127° to the beam at a distance of 32 cm from the target. Absorbers of 8 mm Pb plus 
3 mm Cu and of 3 mm Pb plus 3 mm Cu were placed in front of the detectors at 90° and 
127°, respectively, in the measurements at 7.9 and 9.0 MeV whereas in the measurement at 
13.2 MeV absorbers of 13 mm Pb plus 3 mm Cu and of 8 mm Pb plus 3 mm Cu were used 
for the detectors at 90° and 127°, respectively. A detailed description of the bremsstrahlung 
facility is given in Ref. \d\. 

Spectra of scattered photons were measured for 56 h, 65 h, and 97 h in the experiments 
at 7.9, 9.0, and 13.2 MeV electron energy, respectively. Parts of a spectrum including events 
measured with the two detectors placed at 127° relative to the beam at an electron energy 
of 13.2 MeV are shown in Fig. [H 

A. The photon-scattering method 

In photon-scattering experiments the energy-integrated scattering cross section Ig of an 
excited state at the energy can be deduced from the measured intensity of the respective 
transition to the ground state (elastic scattering). It can be determined relative to the 
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known integrated scattering cross sections Is{E^) 



of states in ^^B 14 : 
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Here, I^{E^,9) and I^{E^,9) denote the measured intensities of a considered ground- 
state transition at E^ and of a ground-state transition in ^^B at E^, respectively, observed 
at an angle 9 to the beam. W{E^,9) and W{E^,9) describe the angular correlations of 
these transitions. The quantities A^jv and are the numbers of nuclei in the ^°Zr and ^^B 
targets, respectively. The quantities ^-y^E^) and stand for the photon fluxes at the 

energy of the considered level and at the energy of a level in ^^B, respectively. 

The integrated scattering cross section is related to the partial width of the ground-state 
transition Fq according to 



where cr^-y is the elastic scattering cross section, E^, Jx and T denote energy, spin and total 
width of the excited level, respectively, and Jq is the spin of the ground state. 

For the determination of the level widths one is faced with two problems. First, a consid- 
ered level can be fed by transitions from higher-lying states. The measured intensity of the 
ground-state transition is in this case higher than the one resulting from a direct excitation 
only. As a consequence, the integrated cross section deduced from this intensity contains 
a part originating from feeding in addition to the true integrated scattering cross section: 
Ig+f = Is + If- Furthermore, a considered level can deexcite not only to the ground state, 
but also to low- lying excited states (inelastic scattering). In this case, not all observed 7 
transitions are ground-state transitions. To deduce the partial width of a ground-state tran- 
sition Fq and the integrated absorption cross section one has to know the branching ratio 
^0 = Tq/F. If this branching ratio cannot be determined, only the quantity Fq/F can be 
deduced (cf. Eq. (E])). 

Spins of excited states can be deduced by comparing experimental ratios of intensities 
measured at two angles with theoretical predictions. The optimum combination are angles 
of 90° and 127° because the respective ratios for the spin sequences — 1 — and 0-2-0 
differ most at these angles. The expected values are iy(90°)/iy(127°)o-i-o = 0.74 and 
IV(90°)/iy(127°)o-2-o = 2.18 taking into account opening angles of 16° and 14° of the 
detectors at 90° and 127°, respectively. 

B. Detector response and photon flux 

For the determination of the integrated scattering cross sections according to Eq. ([T]) 
the relative efficiencies of the detectors and the relative photon flux are needed. For the 
determination of the dipole-strength distribution described in Sec. IIIII the experimental 
spectrum has to be corrected for detector response, for the absolute efficiency, and for 
the absolute photon flux, for background radiation, and for atomic processes induced by 
the impinging photons in the target material. The detector response has been simulated 
using the program package GEANT3 The reliability of the simulation was tested by 
comparing simulated spectra with measured ones as described in Refs. 0, 
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FIG. 2: Absolute photon flux at the target deduced from intensities of transitions in ^^B (circles) 
using the calculated efficiency and the relative photon flux calculated as described in the text (solid 
line) . 



The absolute efficiencies of the HPGe detectors were determined experimentally up to 
1333 keV from measurements with ^^Na, ^''Co, ^^'^Ba, and ^^"^Cs calibration sources. The 
efficiency curve calculated with GEANT3 was scaled to ffi the absolute experimental values. 
To check the simulated efficiency curve at higher energies we used an uncalibrated ^^Co 
source produced in-house using the ^^Fe(p, n)^^Co reaction. The relative efficiency values 
obtained from the ^^Co source adjusted to the calculated value at 1238 keV are consistent 
with the simulated curve as is illustrated in Ref. 0]. 

The absolute photon flux was determined from intensities and known integrated scattering 
cross sections of transitions in ^^B which was combined with the ^°Zr target (cf. Sec. [Tll)- For 
interpolation, the photon flux was calculated using a code 1^ based on the approximation 
given in Ref. [13] and including a screening correction according to Ref. [l8|. In addition, 
the flux was corrected for the attenuation by the beam hardener. This flux was adjusted 
to the experimental values as is shown in Fig. [21 The transition from the level at 7285 keV 
could not be used because of superposition by another transition. 



C. Experiments at various electron energies 

Measurements at various electron energies allow us to estimate the influence of feeding on 
the integrated cross sections. Ratios of the quantities Ig+f obtained for levels in ^°Zr from 
measurements at different electron energies are shown in Fig. |3l The plotted ratios reveal 
that (i) only levels below ^ Q MeV are influenced considerably by feeding and (ii) levels 
in the range of i^^^ ~ 4 to 6 MeV are mainly fed by levels above ^ 9 MeV. Transitions 
found in the measurement at i?g'° = 7.9 MeV are assumed to be ground-state transitions. 
Transitions additionally observed up to 7.9 MeV in the measurements at 9.0 MeV and 13.2 
MeV are therefore considered as inelastic transitions from high-lying to low-lying excited 
states. In the same way, transitions in the range from 7-9 MeV observed in the measurement 
at 9.0 MeV are considered as ground-state transitions, whereas transitions in this energy 
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range found additionally at 13.2 MeV are considered as branchings to low-lying states. By 
comparing the measurements in this way, transitions from high-lying to excited low-lying 
states may be filtered out. The remaining transitions, assumed as ground-state transitions, 
have been used to derive the corresponding level energies which are listed in Table [I] together 
with spin assignments deduced from angular distributions of the ground-state transitions, 
integrated scattering cross sections, and the quantities Fq/F. 

TABLE I: Levels assigned to ^°Zr. 
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41(9) 


341 (73) [8] 


9843.4(6) 


0.77(16) 1 


84(15) 


702(126) 


9855.5(8) 


0.55(17) 1 


58(13) 


488(109) 


9872.4(4) 


0.49(12) 1 


126(24) 


1067(207) 


9890.7(13) 


0.7(5) (1) 


81(34) 


686(289) 


9901.9(13) 




70(28) 


592(242) [8] 


9932.1(12) 


0.43(14) 1 


123(39) 


1053(330) 


9962.8(5) 


0.57(14) 1 


172(37) 


1479(322) 


9984.1(11) 




69(34) 


593(297) [8] 


10004.2(10) 


0.42(19) 1 


70(16) 


606(137) 


10019.6(11) 


0.65(16) 1 


94(16) 


819(138) 


10031(2) 




69(16) 


599(136) [8] 
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10042.9(4) 


0.70(8) (r 


-)[7] 


316(36) 


2762(310) 


10083.8(6) 


0.78(13) 


1 


93(13) 


823(118) 


10094.2(7) 


0.74(19) 


1 


83(14) 


732(120) 


10104.9(12) 


0.7(3) 


(1) 


49(13) 


433(118) 


10123.7(18) 


0.5(3) 


1 


137(100) 


1219(886) 


10146.8(9) 


0.57(27) 


1 


46(13) 


409(114) 


10163.4(8) 


0.50(23) 


1 


60(16) 


540(147) 


10193.0(5) 


0.60(12) 


1 


168(25) 


1514(220) 


10216.8(10) 


0.57(24) 


1 


76(18) 


693(159) 


10233(4) 






47(38) 


427(348) [8] 


10241(2) 


0.9(3) 


(1) 


79(34) 


723(305) 


10260.9(11) 






23(6) 


212(53) [8] 


10270.0(7) 






34(9) 


308(81) [8] 


10286.2(6) 


0.64(19) 


1 


42(9) 


388(79) 


10298.3(10) 


0.8(3) 


(1) 


32(8) 


290(69) 


10306.6(9) 


0.71(21) 


1 


46(9) 


426(80) 


10315.1(4) 


0.63(12) 


1 


103(14) 


952(126) 


10334.9(6) 


0.68(21) 


1 


51(10) 


470(93) 


10361(2) 


0.86(29) 


(1) 


54(14) 


504(129) 


10376.8(4) 


0.46(7) 


1 


240(28) 


2240(259) 


10402.5(9) 


0.72(14) 


1 


85(16) 


799(150) 


10494.5(11) 


0.7(3) 


(1) 


43(10) 


411(92) 


10507.9(8) 


0.6(3) 


1 


49(10) 


467(96) 


10524.6(4) 


0.72(10) 


1 


143(18) 


1376(175) 


10595.0(7) 


0.80(19) 


1 


92(14) 


901(136) 


10618.7(8) 


0.74(23) 


1 


67(12) 


651(119) 


10638.5(9) 


0.41(20) 


1 


59(12) 


575(122) 


10682.2(6) 


0.54(16) 


1 


42(10) 


417(95) 


10713.2(12) 


0.8(4) 


(1) 


37(20) 


369(195) 


10728.2(11) 


0.67(20) 


1 


102(32) 


1022(316) 


10827.1(5) 


0.82(15) 


1 


105(16) 


1068(167) 


10914(2) 


0.79(21) 


(1) 


113(21) 


1168(214) 


10957(2) 


0.35(8) 


1 


118(19) 


1224(202) 


10987.0(10) 


0.66(13) 


1 


161(23) 


1691(242) 


11044(2) 






49(17) 


522(179) [8] 


11094.2(15) 






70(10) 


743(106) [8] 


11108.0(16) 






39(8) 


422(84) [8] 


11120.4(9) 


0.60(16) 


1 


92(17) 


992(179) 


11129.2(17) 






57(18) 


611(196)[8] 


11140(2) 






57(9) 


612(95) [8] 


11232.4(7) 


0.45(14) 


1 


88(13) 


963(146) 


11243.2(6) 


0.58(15) 


1 


92(14) 


1010(152) 


11337.7(6) 


0.85(13) 


1 


91(15) 


1010(167) 


11417.5(7) 


0.8(4) 


(1) 


108(25) 


1217(287) 


11452.2(10) 


0.42(12) 


1 


132(28) 


1501(316) 


11479.7(8) 


0.58(15) 


1 


191(33) 


2181(374) 
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11501(3) 






66(37) 


756(425) [8] 


11510(7) 






33(15) 


382(172) [8] 


11531(2) 


0.42(25) 


1 


74(30) 


848(346) 


11627.9(9) 






44(16) 


518(182)[8] 


11651.5(8) 


1.0(3) 


(1) 


48(16) 


564(185) 


11777.4(10) 


0.46(22) 


1 


124(40) 


1489(479) 


11788(3) 


0.44(29) 


1 


73(36) 


885(440) 


11963.3(18) 


0.8(3) 


(1) 


68(14) 


845(179) 


11984(2) 


0.57(18) 


1 


57(13) 


715(167) 


12020.6(8) 


0.67(14) 


1 


155(21) 


1942(260) 


12067.8(9) 


0.63(17) 


1 


124(19) 


1570(237) 


12208.3(12) 


0.40(15) 


1 


72(16) 


929(214) 


12243.6(14) 


0.57(19) 


1 


62(15) 


799(188) 


12496.3(18) 






87(18) 


1181(244)[8] 


12880.3(10) 






11(3) 


164(46) [8] 



In general, the complete level scheme may be constructed by searching for combinations 
of 7-ray energies that fit another 7-ray according to the Ritz principle. However, if one 
includes all observed transitions one finds a huge number of possible two-fold or even three- 
fold combinations Because we cannot derive complete information about the many 
possible branching transitions, we are not able to determine the branching ratios of the 
ground-state transitions from the present data. Therefore, we applied statistical methods 
to estimate the intensities of branching transitions and thus to correct the dipole-strength 
distribution. 



[1] Excitation energy. The uncertainty in parentheses is given in units of the last digit. This energy was 
deduced from the 7-ray energy measured at 127° to the beam by including a recoil and Doppler correction. 

[2] Ratio of the intensities of the ground-state transitions measured at angles of 90° and 127°. The expected 
values for an elastic pure-dipole transition (spin sequence — 1 — 0) and for an elastic quadrupole transition 

(spin sequence — 2 — 0) are 0.74 and 2.18, respectively. 
[3] Spin and parity of the state. 

[4] Ratio of integrated scattering-Kfeeding cross sections deduced at different electron energies. The deviation 

from unity is a measure of feeding. 
[5] Integrated scattering cross section. The values up to the level at 6875 keV were deduced from the 

measurement at E'g'" = 7.9 MeV, those from the level at 6960 keV up to the level at 8832 keV were 
deduced from the measurement at E^™ = 9.0 MeV, and the values given for levels at higher energies were 
deduced from the measurement at E^"^ — 13.2 MeV. 
[6] Partial width of the ground-state transition Fq multiplied with the branching ratio bo = Tq/T. The values 
deduced from the given 1^ for states up to the 5504 keV state exceed the values of previous studies [l^] 
outside their uncertainties which may indicate that the present values are influenced by feeding. These 
values are therefore not given in this table. An estimate of the partial width Tq can be obtained from Fq/F 
after correction for a mean branching ratio &o = To/F = 0.9(1)% for resolved ground-state transitions (cf. 
SecHnj. 

[7] Spin and parity have been known from previous work (cf. Sec. 

[8] Value deduced under the assumption of J^r = 1. 

[9] This transition is superimposed by a transition in ^^B. 
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III. DETERMINATION OF THE DIPOLE-STRENGTH DISTRIBUTION 



In the following we describe the statistical methods used to estimate the intensity distri- 
bution of branching transitions and the branching ratios of the ground-state transitions. 

By using simulations as described in Sec. IIIBI we corrected the spectrum including the 
two detectors at 127°, measured during the irradiation of the ^°Zr target at = 13.2 
MeV. In a first step spectra of the ambient background adjusted to the intensities of the 
1460.5 keV transition (decay of ^°K) and 2614.9 keV transition (decay of ^os^jj-j ^^le in- 
beam spectrum were subtracted from the measured spectrum. It turned out that transitions 
following (n, 7) reactions in the HPGe detectors and in surrounding materials are negligibly 
small and thus, did not require correction. To correct the spectrum for detector response, 
spectra of monoenergetic 7 rays were calculated in steps of 10 keV by using GEANT3. 
Starting from the high-energy end of the experimental spectrum, the simulated spectra were 
subtracted sequentially. The resulting spectrum including the two detectors at 127° is shown 
in Fig. m 

The background produced by atomic processes in the ^^Zr target was obtained from a 
GEANT3 simulation using the absolute photon flux deduced from the intensities of the 
transitions in ^^B (cf. Fig. [2]). The corresponding background spectrum multiplied with 
the efficiency curve and with the measuring time is also depicted in Fig. HI As can be seen 
in Fig. m the continuum in the spectrum of photons scattered from ^°Zr is clearly higher 
than the background by atomic scattering. This continuum is formed by a large number of 
non-resolvable transitions with small intensities which are a consequence of the increasing 
nuclear level density at high energy and of Porter- Thomas fluctuations of the decay widths 
TqI in connection with the finite detector resolution (e.g. AE f» 7 keV at E"^ ~ 9 MeV). 



^7 

90'- 



The relevant intensity of the photons resonantly scattered from Zr is obtained from a 
subtraction of the atomic background from the response-corrected experimental spectrum. 
The remaining intensity distribution includes the intensity contained in the resolved peaks as 
well as the intensity of the "nuclear" continuum. The scattering cross sections a^^' derived 
from this intensity distribution for energy bins of 0.2 MeV from the full intensity distribution 
are shown in Fig. [51 These values are compared with those given in Table [I] for resolved 
transitions in ^°Zr. One sees that the two curves have similar structures caused by the 
prominent peaks. However, the curve including also the continuum part of the spectrum 
contains altogether a strength that is by a factor of about 2.7 greater than the strength of 
the resolved peaks only. The full intensity distribution (resolved peaks and continuum) and 
the corresponding scattering cross sections shown in Fig. [5] contain ground-state transitions 
and, in addition, branching transitions to lower-lying excited states (inelastic transitions) as 
well as transitions from those states to the ground state (cascade transitions). The different 
types of transitions cannot be clearly distinguished. However, for the determination of the 
photoabsorption cross section and the partial widths Tq the intensities of the ground-state 
transitions are needed. Therefore, contributions of inelastic and cascade transitions have 
to be subtracted from the spectra. We corrected the intensity distributions by simulating 
7-ray cascades 2^ from the levels with J = 0, 1, 2 in the whole energy range analogously 



to the strategy of the Monte-Carlo code DICEBOX [2l||. In these simulations, 1000 nuclear 



realizations starting from the ground state were created with level densities derived from 



experiments [22|. We applied the statistical methods also for the low-energy part of the 
level scheme instead of using experimentally known low-lying levels in ^°Zr because this 
would require the knowledge of the partial decay widths of all transitions populating these 
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fixed levels. Fluctuations of the nearest-neighbor-spacings were taken into account according 



to the Wigner distribution (see, e.g., Ref. [23|). The partial widths of the transitions to 



low-lying levels were assigned using a priori known strength functions for El, Ml, and E2 
transitions. Fluctuations of the partial widths were treated the applying Porter-Thomas 
distribution 1^ • 



In the calculations, the recently published parameters for the Back-Shifted Fermi-Gas 



(BSFG) model obtained from fits to experimental level densities 22!], cl = 8.95(41) MeV ^ 



and El = 1.97(30) MeV, were used. In the individual nuclear realizations, the values of a 
and El were varied within their uncertainties. As usual in the BSFG model, we assumed 
equal level densities for states with positive and negative parities of the same spin j22| . This 
assumption has recently been justified by the good agreement of level densities predicted 
by the BSFG model with experimental level densities of 1+ states in the energy range 
from 5 to 10 MeV obtained from the ^''Zr(^IIe,t)^°Nb reaction [13] and with experimental 
level densities of 2+ and 2~ states in ^°Zr studied in the ^°Zr(e, e') and ^°Zr(p,p') reactions 



25l |. The extended analysis of the ^°Zr(^He,t)^°Nb reaction in Ref. [25| indicates however 
fluctuations of the level density of 1"^ states in ^°Nb around the predictions of the BSFG 
model. 



For the El, Ml, and E2 photon strength functions Lorentz parametrizations [26 



were 



used. The parameters of the Lorentz curve for the El strength were taken from a fit to 



(7, n) data [27|] and are consistent with the Thomas-Reiche-Kuhn sum rule (TRK) fcroF 



60NZ/A MeV mb 2^. The parameters for the Ml and E2 strengths were taken from global 
parametrizations of Ml spin-flip resonances and E2 isoscalar resonances, respectively [2^ . 

Spectra of 7-ray cascades were generated for groups of levels in 100 keV bins in each 
of the 1000 nuclear realizations. For illustration, the distributions resulting from 10 indi- 
vidual nuclear realizations populating levels in a 100 keV bin around 9 MeV are shown in 
Fig. O which reflect the influence of fluctuations of level energies and level widths. Because 
in the nuclear realizations the levels were created randomly starting from the ground state 
instead of starting with the known first excited state at 2.2 MeV, the distribution of the 
branching transitions continues to the energy bin of the ground-state transitions. These 
spectra resemble qualitatively the ones measured in an experiment on ^°Zr using tagged 
photons Starting from the high-energy end of the experimental spectrum, which con- 
tains ground-state transitions only, the simulated intensities of the ground-state transitions 
were normalized to the experimental ones in the considered bin and the intensity distribu- 
tion of the branching transitions was subtracted from the experimental spectrum. Applying 
this procedure step-by-step for each energy bin moving toward the low-energy end of the 
spectrum one obtains the intensity distribution of the ground-state transitions. Simultane- 
ously, the branching ratios of the ground-state transitions are deduced for each energy 
bin A. In an individual nuclear realization, the branching ratio is calculated as the ra- 
tio of the sum of the intensities of the ground-state transitions from all levels in A to the 
total intensity of all transitions depopulating those levels to any low-lying levels including 
the ground state 0, [1]. By dividing the summed intensities in a bin of the experimental 
intensity distribution of the ground-state transitions by the corresponding branching ratio 
we obtain the absorption cross section for a bin as = a^^/b^. Finally, the absorption 
cross sections of each bin were obtained by averaging over the values of the 1000 nuclear 
realizations. For the uncertainty of the absorption cross section a la deviation from the 
mean has been taken. 

To get an impression about the branching ratios, the individual values of 10 nuclear 
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realizations are shown in Fig. [71 The mean branching ratio of the 1000 reahzatios decreases 
from about 80% for low-lying states, where only few possibilities for transitions to lower- 
lying states exist, to about 65% at the neutron-separation energy. Toward low energy the 
uncertainty of increases due to level-spacing fluctuations and the decreasing level density. 
The large fluctuations below about 6 MeV make these values useless. Note that the mean 
branching ratio is not representative for transitions with large intensities like the resolved 
transitions given in Table [H It turns out from the simulations that the branching ratios of 
transitions with partial widths Fq like the ones given in Table [T] are in the order of bo ~ 85% 
to 99% (cf. also Ref. Q). 

The photoabsorption cross sections obtained for ^°Zr are compared with the data of a 

Sevious experiment with monoenergetic photons in the energy range from 8.4 to 12.5 MeV 
in Fig. [HI The data of Ref. [1] contain corrections in form of additional constant partial 
widths for branching and for pro ton emission, i.e. for the (7,p) channel. Cross sections 
obtained in (7, n) experiments |27j and cross sections calculated for the (7, p) reaction using 
the Talys code [30(] with the Lorentz model for the El strength function [26| are also shown 
in Fig. E 

To test the influence of the El strength function on the results of the simulation and to 
check the consistency between input strength function and deduced photoabsorption cross 
section we performed cascade simulations with modified El strength functions. In one case, 
we applied a Lorentz curve with an energy- dependent width T{E) ~ In another case 
we assumed a FDR on top of the Lorentz curve with constant width. The FDR had a 
Breit-Wigner shape with its maximum at 9 MeV and a width of 2 MeV. The cross section 
of the FDR was chosen 2.2% of the TRK (see discussion in Sec. lIVBp . Cross sections 
corresponding to these input strength functions and the obtained absorption cross sections 
are shown in Fig. [91 One sees that the resulting absorption cross sections differ only slightly 
from each other and also from the ones obtained by using the Lorentz curve with constant 
width (cf. Fig. [8]) and overlap within their uncertainties. This shows that the resulting cross 
sections are not very sensitive to modifications of the input strength functions. However, the 
comparison of the resulting cross sections with the input strength functions shows that there 
is consistency in the case of the Lorentz curve with a FDR in addition, whereas there is a large 
discrepancy in the case of the Lorentz curve with energy-dependent width. Summarizing, 
the Lorentz curve with a constant width and this curve with a FDR in addition deliver 
consistency between input and output strength functions and similar results. 

The total photoabsorption cross section has been deduced by combining t he p resent 
(7,7') data with the (7,77,) data of Ref. 27| and calculated (7,p) cross sections [30|. This 
total cross section is compared with the Lorentz curve described above, which is currently 
used as a standard strength function for the calculation of reaction data [2^ , in Fig. [TOl It 
can be seen that the experimental cross section includes extra strength with respect to the 
approximation of the low-energy tail of the GDR by a Lorentz curve in the energy range 
from 6 to 11 MeV and is in better agreement with a strength function including a FDR on 
top of the Lorentz curve (cf. Fig. To investigate the nature of this extra strength we 
have performed model calculations which are discussed in the following section. 



13 



IV. QUASIPARTICLE-PHONON-MODEL CALCULATIONS FOR THE iV = 50 
ISOTONES 88Sr AND ^^Zr 

The experimental data on the dipole response of ^°Zr discussed in the present work reveal 
a strong enhancement of the photoabsorption cross section (see Fig. [TU]) in the energy range 
i?^ = 6 - 11 MeV and a considerable deviation of its shape from the Lorentz-like strength 
function used to adjust the GDR. Similar dipole resonance structures have been detected 
also in our recent study of ^^Sr 0]. In the following we investigate the nature of the dipole 
strength in the iV = 50 isotones ^^Sr and ^°Zr in the framework of the extended quasiparticle- 
phonon model (QPM) jsi, 3^. The approach is based on a Hartree-Fock-Bogoljubov (HFB) 



description of the ground state [31| by using a phenomenological energy-density functional 
(EDF) {32! . The excited states are calculated within the QPM j33| . 



A. The QPM model 

The model hamiltonian is given by: 



H = HMF + Hf,+Hll, + HZ (3) 



and resembles in structure that of the traditional QPM model |33| but in detail differs in 
the physical content in important aspects as discussed in Ref. 31, 33] ■ The HFB term 
Hmf = Hsp + Hpair contains two parts. Hgp describes the motion of protons and neutrons 
in a static, spherically-symmetric mean field and Hpair accounts for the monopole pairing 
between isospin-identical particles with phenomenologically adjusted coupling constants, 



fitted to data [3J]. The mean- field hamiltonian compromises microscopic HFB effects |32l. |35 
and phenomenological aspects like experimental separation energies and, when available, also 
charge and mass radii. For numerical convenience the mean-field potential is parametrized 
in terms of a superpositions of Wood-Saxon (WS) potentials with adjustable parameters 

The remaining three terms present the residual interaction Hres = H^m + -^5m + -^m- 
As typical for the QPM, we use separable multipole-multipole and spin-multipole -^5^ 
interactions both of isoscalar and isovector type in the particle-hole channel and a multipole- 
multipole pairing in the particle-particle channel. The isoscalar and isovector coupling 
constants are obtained by a fit to energies and transition strengths of low-lying vibrational 
states and high-lying GDRs j32|, 3^ . 



The excited nuclear states are described by quasiparticle-random-phase-approximation 
(QRPA) phonons. They are defined as a linear combination of two-quasiparticle creation 
and annihilation operators: 

= In - ^';'^A,(jj')) , (4) 

jj' 

where j = {nljrriT) is a single-particle proton or neutron state, A'l^ and Ax^ are time-forward 
and time-backward operators, coupling proton and neutron two-quasiparticle creation or 
annihilation operators to a total angular momentum A with projection n by means of the 
Clebsch-Gordan coefficients C^^^-,^, = {jmj'm'\Xfj,). The time reversed operator is defined 

as Ax^ = (-)^-M, 
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We take into account the intrinsic fermionic structure of the phonons by the commutation 
relations: 



[Qx/ii, Qx'u'i'] 



3J 



mj2m2 



3] 32 
mm' 1712 



_('_NA+A'+^+/i' Ai X'i' ^X-fi f-iX'-n' 
\ J ^332^3'h 3fnj2m2^j'm'j2m2' 

The QRPA phonons obey the equation of motion 



(5) 



(6) 



which solves the eigenvalue problem, i.e. it gives the excitation energies Ei and the time- 
forward and time-backward amplitudes [sH il^j^j^ and ^Pjlj^, respectively. 

In terms of the QRPA phonons the model hamiltonian can be expressed as 



H - Hph + Hqph - ^ ^XiQt^iQx^ii 
X/j-i 



(7) 



+ 2 XI C'AiVA2M2^Aiii(^3«3) 

AiA2A3jii2i3MiA'2M3 



[^AiMljl^A2M2i2^A3-M3j3 + ^■'^■]- 

The first part in the Eq. (I7j) refers to the harmonic part of nuclear vibrations, while the 
second one is accounting for the interaction between quasiparticles and phonons. The latter 
allows us to go beyond the ideal harmonic picture by including effects due to the non-bosonic 
features of the nuclear excitations. 

The model hamiltonian is diagonalized assuming a spherical O"*" ground state which leads 
to an orthonormal set of wave functions with good total angular momentum JM. For even- 
even nuclei we extend the approach to anharmonic effects by considering wave functions 
including a mixture of one-, two-, and three-phonon components 37| : 



+ 

JMi 



JM 



(8) 



Alii 
A2i2 



-^in^2»2 



2-AinA2i2/(j X rn+ . 0+ .1 (g)Q 



A3M3«3 



JM 



where i> labels the number of the excited states. The energies and the phonon-mixing 
amplitudes R, P, and T are determined by solving the extended multi-phonon equations of 
motions. 
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TABLE II: HFB results obtained with the phenomenological EDF for total binding energy per 
nucleon, the measured binding energy [34], proton and neutron root-mean-square radii, and neutron 
skin thickness, Eq. ([9]), for the two = 50 isotones. 





^e.p (MeV) 


^t.e (MeV) 


(fm) 


^/rl (fm) 


8r (fm) 


88Sr 


-8.733 


-8.766 


4.185 


4.291 


0.106 


902r 


-8.710 


-8.692 


4.229 


4.303 


0.074 



As discussed in Refs. 32|, |38|] we obtain important additional information from the analysis 
of the spatial structure of the nuclear response on an external electromagnetic field. That 
information is accessible by considering the one-body transition densities Sp{f), which are 
the non-diagonal elements of the nuclear one-body density matrix. The transition densities 
are obtained by the matrix elements between the ground state |\E'j) = | JjMj) and the excited 
states |\E'j) = \JfMf). We identify in QPM |JjMj) = |0) with phonon vacuum and obtain 
the excited states by means of the QRPA state operator, Eq. (jl]), \ JfMf) = Qx^i\0). The 
analytical procedure for the calculation of the transition densities is presented in details in 
Ref. [si. 



B. Calculation of dipole excitations in the A = 50 isotones 

In Ref. 3^, a close relationship between skins in the ground-state matter-density distri- 
butions and the low-energy dipole response was pointed out. It was shown, that the skin 
thickness, defined by 

Sr = ^< rl > - J<rl> (9) 



is indeed directly connected to the non-energy-weighted dipole sum rule, thus being com- 
plementary to the well established relation of the GDR strength and the energy-weighted 
Thomas- Reiche-Kuhn dipole sum rule. From the calculated neutron and proton ground state 
densities displayed in Fig. [H] it is seen that the two A^ = 50 isotones considered here, ^^Sr 
and ^°Zr, show a neutron skin. The total binding energies and radii are given in Table [Tll 
Of special importance for our investigations are the surface regions, where the formation of 
a skin takes place. For the investigated ®^Sr and ^°Zr nuclei, the neutron skin decreases with 
decreasing N/Z, i.e. when going from ®^Sr to ^°Zr. 

The QRPA calculations in these nuclei result in a sequence of low- lying one-phonon 
dipole states of almost pure neutron structure, located below the particle threshold. The 
proton contribution in the state vectors is less than 5%. These states have been related 
to the oscillations of weakly bound neutrons from the surface region, known as neutron 



PDR (see Refs. [3l|, [32, l38|, |39|, |40| and Refs. therein). The calculations of the dipole 
response in A^ = 50 up to 22 MeV are presented in Fig. O The El transition matrix 
elements were calculated with recoil-corrected effective charges qn = —Z/A for neutrons and 



Qp = N/A for protons [32|, respectively. To compare the predicted dipole strengths with 
the experimental values we calculated the integrated absorption cross sections according to 
fa^dE = 4:.03E^B{E1, 0+ 1") MeV mb with E^ in MeV and B{E1) in e^fm^ |4l|. When 
comparing the QRPA results with experimental cross sections one has the general problem 
of comparing discrete values (cf. Fig. [T^ with gradually changing averages of energy bins 
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(cf. Fig. [T0|) . Therefore, the QRPA values are usually folded with Lorentz curves that can 
be considered as a simulation of the damping of the GDR due to higher-order configurations 
neglected in QRPA. In the present case, we folded the discrete QRPA cross sections with 
Lorentz curves of 3.0 MeV width, which is a typical value for the damping width. The 
results are compared with the experimental cross sections in Fig. [10] for ^°Zr and in Fig. [T3] 
for ^^Sr. The QRPA calculations predict extra strength with respect to the Lorentz curve 
in qualitative agreement with the experimental findings. However, the experimental cross 
sections show a clear separation of the PDR and the GDR regions around 12 MeV, whereas 
the calculations seem to predict an excess of strength up to about 14 MeV. 

The nature of the dipole strength can be examined by analyzing dipole transition densities 
that act as a clear indicator of the PDR and allow us to distinguish between the different 
types of the dipole excitations 3^. The QRPA proton and neutron dipole transition densities 
for ^^Sr and ^°Zr were summed over selected energy regions and are shown in Fig. O They 
can be directly related to the spectrum of 1~ states in Fig. [T21 In the region of -E^ < 9 
MeV we find an in-phase oscillation of protons and neutrons in the nuclear interior, while 
at the surface only neutrons contribute. Because these features are characteristic for dipole 
skin modes, we identify these states as part of the PDR. The states in the region of from 
9 MeV to 9.5 MeV carry a different signature and their contribution to the total proton 
and neutron transition densities in the energy region < 9.5 MeV changes completely the 
behavior of the proton/neutron oscillations related to the PDR. At these higher energies the 
protons and neutrons start to move out of phase being compatible with the low-energy part 
of the GDR. At E^ from 9 to 22 MeV a strong isovector oscillation corresponding to the 
excitation of the GDR is observed. 

In the case of iV = 50 isotones the neutron number is fixed and only the proton number 
changes. This affects the thickness of the neutron skin as well as is seen in Fig. [11] and 
Table HTl Correspondingly, the total PDR strength decreases with increasing proton number. 
According to the criteria discussed above we consider the states in the energy region E^ < 
9 MeV as a part of the PDR. In this energy interval, the QRPA calculations predict a total 
PDR strength of ^5(E1)T= 0.107 e^fm^ in ^oZr and ^5(E1)T= 0.170 e^fm^ in ^^Sr, 
respectively. This accounts for 0.27% and 0.44% of the TRK in ^°Zr and ^^Si, respectively, 
and compares with experimental values of 0.55(2)% and 0.74(2)% derived from resolved 
peaks in ^°Zr and ^^Sr, respectively, whereas the values deduced from the corrected full 
intensities including also the continuum part are remarkably greater and amount to 2.2(4)% 
and 2.4(6)% in ^°Zr and ^'^Sr, respectively. 

Finally, we note that the results are in agreement with the established connection between 
the calculated total PDR strength and the nuclear skin thickness, Eq. ([9]), found in our 
previous investigations in the = 82 isotones and the ^^^Sn isotopes [s^l nuclei. 



In the multi-phonon QPM calculations the structure of the excited states is described by 
wave functions as defined in Eq. ([H]). We now investigate multi-phonon effects using a model 
space with up to three-phonon components, built from a basis of QRPA states with J'" = 
1"*", 1~, 2"*", 3~, 4+, 5~, 6^, 7~, and 8~^. Because the one-phonon configurations up to E^ = 
22 MeV are considered the core polarization contributions to the transitions of the low-lying 
1~ states are taken into account explicitly. Hence, we do not need to introduce dynamical 
effective charges. Our model space includes one-, two- and three-phonon configurations (in 
total 200) with energies up to E^ = 12 MeV. 

The experimental absorption cross sections in ^°Zr in the range from 5 to 11 MeV are 
compared with the corresponding QPM results in Fig. [151 The strength in the region E^ 



17 



= 6-7 McV is related to the PDR as identified from the structure of the involved 1~ 
states. The contribution of the GDR phonons becomes significant at energies > 7.5 MeV 
where a coupling between PDR, GDR, and multiphonon states reflects the properties of the 
dipole excitations at higher energies. At low energy, the calculated values underestimate 
the experimental values derived from the full intensity and are close to the values deduced 
from resolved peaks only, whereas at high energy they are in better agreement with the 
experimental values including the full intensity. The calculated EWSR from 5 to 9 MeV of 
1.3% of the TRK underestimates the experimental value of 2.2(4)%. However, increasing the 
energy range up to 11 MeV the theoretical value of 3.3% is much closer to the experimental 
result of 4.0(3)%, indicating a slightly different distribution of the QPM strength with respect 
to excitation energy. 

V. SUMMARY 

The dipole-strength distribution in ^°Zr up to the neutron-separation energy has been 
studied in photon-scattering experiments at the ELBE accelerator using various electron en- 
ergies. Ground-state transitions have been identified by comparing the transitions observed 
at different electron energies. We identified about 190 levels with a total dipole strength 
of about 180 meV/MeV'^ Spin J = 1 could be deduced from the angular correlations of 
the ground-state transitions for about 140 levels including a dipole strength of about 155 
meV/MeV'\ 

The intensity distribution obtained from the measured spectra after a correction for 
detector response and a subtraction of atomic background in the target contains a continuum 
part in addition to the resolved peaks. It turns out that the dipole strength in the resolved 
peaks amounts to about 27% of the total dipole strength while the continuum contains about 
73%. 

An assignment of inelastic transitions to particular levels and, thus, the determination 
of branching ratios was in general not possible. To get information about the intensities 
of inelastic transitions to low-lying levels we have applied statistical methods. By means 
of simulations of 7-ray cascades intensities of branching transitions coTild be estimated and 
subtracted from the experimental intensity distribution and the intensities of ground-state 
transitions could be corrected in average for their branching ratios. 

A comparison of the photoabsorption cross section obtained in this way from the present 
(7, 7') experiments with (7, n) data shows a smooth connection of the data of the two dif- 
ferent experiments and gives new information about the extension of the dipole-strength 
function toward energies around and below the threshold of the (7, n) reaction. In compari- 
son with a straightforward approximation of the GDR by a Lorentz curve one observes extra 
El strength in the energy range from 6 to 11 MeV which is mainly concentrated in strong 
peaks. 

QPM calculations in *^Sr and ^^Zr predict low-energy dipole strength in the energy region 
from 6 to 10 MeV. The states at about 6 to 7.5 MeV have a special character. Their structure 
is dominated by neutron components and their transition strength is directly related to the 
size of a neutron skin. Their generic character is further confirmed by the shape and structure 
of the related transition densities, showing that these PDR modes are clearly distinguishable 
from the GDR. The dipole states, seen as a rather fragmented ensemble at E^ > 7.5 MeV 
mix strongly with the low-energy tail of the GDR starting to appear in the same region. 
The comphcated structure of these states and the high level densities imposes considerable 
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difficulties for a reliable description of the fragmentation pattern. 

The present analysis shows that standard strength functions currently used for the calcu- 
lation of reaction data do not describe the dipole-strength distribution below the (7, n) 
threshold correctly and should be improved by taking into account the observed extra 
strength. 
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FIG. 3: Ratios of integrated cross sections Ig+f of transitions in ^°Zr obtained at different electron 
energies. 
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FIG. 4: Experimental spectrum of ^^Zr (corrected for room background and detector response) 
and simulated spectrum of atomic background (multiplied with efficiency and measuring time). 




FIG. 5: Scattering cross sections in ^^Zr, derived as a^y = ^^/s/A, not corrected for branching 
and averaged over energy bins of A = 0.2 MeV, as derived from the difference of the experimental 
spectrum and the atomic background shown in Fig. [4] ("peaks + cont.", triangles) and from the 
resolved peaks only ("peaks", circles). 
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FIG. 6: Simulated intensity distribution of transitions depopulating levels in a 100 keV bin around 
9 MeV in ^°Zr. The squares depict the intensities obtained from 10 individual nuclear realizations. 
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FIG. 7: Branching ratios of ground-state transitions as obtained from the simulations of 7-ray 
cascades for ^"^Zr. The squares represent the values of 10 individual nuclear realizations. 
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FIG. 8: (Color online) Photoabsorption cross section deduced from the present (7,7') data for ^''Zr 
after correction for branching transitions (red circles) in comparison with data obtained from an 



experiment with monoenergetic photons (open black circles) [2], (7,?^) data taken from Ref. 27 1 



(green boxes) and with (7,^) cross sections calculated with the Talys code [30|] (blue triangles). 
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FIG. 9: (Color online) Photoabsorption cross sections deduced from the present (7,7') data for 
^°Zr like the red circles in Fig. [51 but with modified El strength functions. The black squares are 
the results for a Lorentz curve with energy-dependent width (black line) and the red circles are 
the results for the assumption of a PDR on top of a Lorentz curve with constant width (pink line). 
The Lorentz curve with a constant width and without PDR is shown as a black dashed line. 
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FIG. 10: (Color online) Total photoabsorption cross section of ^^Zi obtained by combining the 
present (7,7') data, the (7,n) data of Ref. [27| and the (7,^) calculations of Ref. [30(]. The data 
were averaged over 0.5 MeV bins to reduce statistical fluctuations. The black dashed line represents 
a Lorentz curve with the parameters taken from Ref. [27] as also shown in Fig. [9l The blue solid 
line is the result of QRPA calculations discussed in Sec. IIVBI 
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FIG. 12: Calculated dipole one-phonon transition strengths up to Ex = 22 MeV for the isotones 
^^Sr and ^^Zr. 



100 



J3 

E 



10 




5 6 7 8 9 10 11 12 13 14 15 16 17 18 

(MeV) 

FIG. 13: (Color online) Total photoabsorption cross section of ®^Sr as taken from Ref. The 
black dashed line represents a Lorentz distribution with the parameters given in Ref. 0. The blue 
solid line is the result of the present QRPA calculations. 
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FIG. 14: Proton (solid lines) and neutron (dashed lines) transition densities in the = 50 nuclei 
s^Sr and ^^Zr. 
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FIG. 15: (Color online) Experimental cross sections derived from resolved peaks (squares), from 
the full intensity corrected for inelastic transitions and branching ratios (circles) , and cross sections 
predicted by QPM calculations (blue triangles) for ^^Zr, averaged over energy bins of 0.5 MeV. 
Lines are drawn to guide the eye. 
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